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ABSTRACT 

We determine the velocity distribution and space density of a volume complete 
sample of A and F stars, using parallaxes and proper motions from the Hipparcos 
satellite. We use these data to solve for the gravitational potential vertically in the local 
Galactic disc, by comparing the Hipparcos measured space density with predictions 
from various disc models. We derive an estimate of the local dynamical mass density 
of 0.102 ± 0.010 Mq pc" 3 , which may be compared to an estimate of 0.095 M Q pc -3 in 
visible disc matter. Our estimate is found to be in reasonable agreement with other 
estimates by Creze et al. (1998) and Pham (1997), also based on Hipparcos data. We 
conclude that there is no compelling evidence for significant amounts of dark matter 
in the disc. 

Key words: Galaxy: kinematics and dynamics - (Galaxy:) solar neighbourhood - 
Galaxy: structure - dark matter 



CLf 1 INTRODUCTION 



The European Space Agency's Hipparcos satellite has pro- 
duced very accurate distances and proper motions for a com- 
plete set of nearby stars. These data permit a reassessment 
of the amount of matter in the local Galactic disc because 
of two significant improvements. Firstly, the kinematics and 
vertical density distribution of stars used to trace the disc 
potential can be determined with much higher accuracy, per- 
mitting a better determination of the total amount of gravi- 
tating matter in the local disc. Secondly, Hipparcos improves 
the measurement of the local Luminosity Function, so that 
the amount of matter in the disc in visible components can 
be better estimated. Both these improvements lead to a bet- 
ter evaluation of any difference between the amount of vis- 
ible disc matter (in gas and dust, stars, stellar remnants 
and sub-stellar objects) and the total amount of gravitating 
matter. 

Any difference between these quantities would imply 
that there remain mass components in the disc which have 
not yet been directly observed, i.e. disc dark matter. The 
first proposal that there might be a lot of such unobserved 
matter dates to Oort (1932, 1960). Oort used the now clas- 
sical method of solving the combined Poisson-Boltzmann 
equations for the kinematical and density distribution of a 
population of stars, assumed to be stationary in the total 
matter distribution of the disc. Oort found that approx- 
imately one third of the local disc mass remained unac- 
counted for. Modern studies can arguably be dated from the 
work of Bahcall (1984a, b,c) who introduced a new method 



of describing the visible disc matter as a series of isother- 
mal components, in a reanalysis and corroboration of Oort's 
results. 

Studies since 1984 fall into two types, with mildly con- 
flicting results. Firstly, the measurement of the kinematics 
and vertical density falloff in the disc of a stellar tracer (with 
the important improvement that both are determined from 
a single sample) from which the local volume density and/or 
local column density of matter in the disc may be derived. 
Kuijken & Gilmore (1989a,b,c), Kuijken (1991) found lit- 
tle evidence for disc dark matter using faint K dwarfs at 
the South Galactic Pole (SGP), whereas Bahcall, Flynn & 
Gould (1992), using bright K giants at the SGP, confirmed 
BahcalPs (1984) claim of dynamically significant disc dark 
matter (over 50 per cent by mass) although with weakened 
statistical significance. The inclusion of the kinematics and 
local density of nearby K giants by Flynn & Fuchs (1994) to 
the sample of Bahcall et al. (1992) reduced the discrepancy 
between visible and total matter; they derived best fits with 
only 20 per cent of the disc mass in dark form. Secondly, 
general star counts in combination with mass models of the 
Galaxy may be used to infer the amount of disc dark matter. 
Bienayme, Robin & Creze (1987), in a study of star counts to 
faint magnitudes, found that the data could be well fitted 
by considering the mass distribution of the known matter 
only. 

Very recently Creze et al. (1998) have used the results 
from the Hipparcos mission to measure the kinematics and 
vertical falloff of complete samples of nearby A-F stars. The 
determination of the local density of gravitating matter then 
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proceeds via the Poisson-Boltzmann equations for the den- 
sity of a tracer component in a self-gravitating disc. Only 
for A and F type stars does the density change sufficiently 
with vertical distance z from the disc mid-plane that the to- 
tal local density of matter can be measured. The very great 
improvement comes because the distances to the tracer stars 
are measured directly for a complete sample via parallaxes, 
rather than via photometric methods as has always been 
the case before Hipparcos. They find that no dark matter 
is required in the local disc, measuring the total amount of 
gravitating matter as po = 0.076 ± 0.015 M Q pc while the 
amount of visible matter is estimated to lie in the range 
Pvis = [0.06,0.10] M Q pc" 3 . 

Another study using Hipparcos data was made by Pham 
(1997), who used a different method to study F stars, where 
he determined scale heights and velocity dispersions for his 
chosen sample and then, assuming isothermality, used the 
2irGh^po = al relation. He finds the local mass density to 
be po = 0.11 ±0.01 M Q pc" 3 . 

In this paper we have also used A-F stars in the Hip- 
parcos catalogue to measure the local density of matter po, 
using an analysis which differs somewhat from that of Creze 
et al. (1998) and Pham (1997). We confirm the Creze et 
al. (1998) and Pham (1997) result of little significant disc 
dark matter. We use the method of von Hoerner (1960) as 
described by Flynn & Fuchs (1994), in which the velocity 
distribution of a tracer population at the mid-plane may be 
integrated in a model of the local disc potential to yield its 
density falloff in the vertical direction. The predicted falloff 
may then be compared with the Hipparcos measured falloff, 
and the disc model evaluated. As is usual, our disc mod- 
els describe the density stratification of known and putative 
dark components in the disc. 

In section 2 we describe the selection of the sample from 
Hipparcos. In section 3 we detail our models of the disc 
mass stratification, and in section 4 we fit the data with 
and without dark matter in the models, finding it unneces- 
sary to invoke disc dark matter. In section 5 we summarize 
and conclude. 



2 THE HIPPARCOS CATALOGUE AND 
SELECTION OF SAMPLE 

In order to measure the local density of matter, we require 
the vertical kinematics (i.e. wo velocity distribution at the 
Galactic midplane, / (|ioo|)) and the vertical density law v(z) 
of a suitable tracer population. Hipparcos parallaxes may be 
used to estimate distances of tracer stars up to 200 pc. Only 
A and F stars develop a useful change in the density as a 
function of height z above the plane within 200 pc, so it is 
to these stars we confine our study. 

The stars are drawn from the Hipparcos Survey (here- 
after the "Survey"), a predefined part of the Hipparcos Cat- 
alogue intended to be a complete, magnitude limited stellar 
sample. For stars of spectral type earlier than or equal to 
G5 this magnitude limit is V < 7.9 + 1.1 sin \b\. 

We select the A-F stars as follows. We select Survey 
stars with —0.2 < B — V < 0.6 and apply an absolute mag- 
nitude cut of 0.0 < Mv < 2.5. These cuts are illustrated 
in Fig. 1. This results in 14342 stars. Proper motions and 
parallaxes are available for all these stars. We further divide 




Figure 1. Hcrtzsprung- Russell diagram of Hipparcos survey stars 
with a-^/ir < 15 per cent. The sample used in this study is con- 
fined within the box. The dashed line divide the stars into the 
"A"-star (above the line) and "F"-star sample. 



this sample into two at an absolute magnitude My — 1.0. 
We refer to these two samples as A and F stars. The first 
sample roughly takes in B5 to A5 stars, while the second 
covers AO to F5 as a result of the rather large dispersion in 
the absolute magnitude to spectral type relation shown by 
Hipparcos (Houk et al. 1997; Jaschek & Gomez 1998). 

The vertical velocity dispersion of stars changes quite 
quickly as a function of absolute magnitude along the upper 
main sequence, as a consequence of the age-velocity relation 
(Wielen 1977). For this reason, dividing the sample into A 
and F makes the analysis of the velocity-density information 
in Hipparcos considerably simpler. Furthermore, since the A 
star sample is likely to be considerably younger than the F 
star sample, there is a possibility that the stars are so young 
that they are not yet fully mixed into the Galactic potential, 
a key assumption in applying the method. By working with 
an A and an F star sample we have some chance of detecting 
such effects. Ideally, one would prefer to be using old stars, 
such as dwarfs below the turnoff or K giants. We cannot use 
the dwarfs since they are well below the faintness limit of 
the Survey part of Hipparcos. The K giants develop only a 
small change in density with vertical distance within 200 pc. 

The A star sample contains 4448 stars. We remove stars 
further than 200 pc, leaving 2026 stars. The F star sample 
contains 9894 stars. These stars are intrinsically fainter than 
the A stars, and we must apply a 100 pc distance limit in 
order to ensure that the sample remains complete. Applying 
this cut leaves 3080 stars. 

For these two samples we use the Hipparcos proper mo- 
tions and parallaxes to determine the distribution of vertical 
velocity and change in density with height, described in the 
next two subsections. 
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2.2 Vertical velocity distribution 
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Figure 2. Lcfthand panels show the number of stars in 10 pc 
slices in the samples before (solid) and after (dotted) corrections 
for distance measurement errors and reddening for A-stars (upper 
panel) and F-stars (lower panel). The righthand panels show the 
relative vertical distance error as a function of distance to the 
sun, z — 20, for the two samples. 



2.1 Stellar density laws 

The density distribution of the tracer samples as a function 
of vertical distance from the galactic plane is determined 
within a cylinder centered on the sun. The radius of the 
cylinder is 100 pc for the F-star sample, and 200 pc for the 
A-star sample. Galactic (x, y, z) coordinates are determined 
for each star from the parallax and position on the sky, and 
the number of stars in 10 pc slices in z is determined. These 
samples are sufficiently local that the effects of dust absorp- 
tion are very minor. We correct for extinction effects using 
the extinction model of Hakkila et al. (1997), which is a 
synthesis of a number of published studies. 

The number density falloff and mean distance error as 
a function of vertical distance from the Sun \z — z® \ is shown 
for the A-stars in the upper panels of Fig. 2 and for the F- 
stars in the lower panels. For most of the bins, the error in 
the distances to the stars is much less than the bin width, 
meaning that corrections to the observed density are small. 
For the A-stars the reddening correction has a larger im- 
pact than the distance errors, but only results in a change 
of the normalization of the counts and does not effect the 
slope, which is what we are concerned with here. To illus- 
trate this, the solid histograms in the left panels show the 
observed density law in \z — zq\, while the dotted lines show 
the extinction corrected density law. The corrections have 
been calculated using Monte-Carlo methods, in which large 
numbers of stars are simulated on the sky using the galactic 
model of Holmberg, Flynn & Lindegren (1997). We simu- 
late the observation of stars by Hipparcos in the models, 
including the errors as a function of apparent magnitude 
and position on the sky, and this allows us to calculate the 
small correction between the observed and true density. 



Hipparcos gives us for the first time large numbers of ac- 
curate parallaxes and proper motions for a kinematically 
unbiased sample. Unfortunately, the original plan to obtain 
complementary radial velocities measurements was not im- 
plemented although the situation will improve in the near 
future (Udry et al. 1997). Full space velocities (u,v,w) are 
unobtainable for most of the stars, and any subsample of 
stars that have radial velocities would be kinematically bi- 
ased towards high velocity stars. To get an unbiased esti- 
mate of the vertical velocities w we sacrifice the radial ve- 
locity information and work exclusively with the Hipparcos 
parallaxes and proper motions. 

We use the standard galactic triad with principal direc- 
tions x towards the Galactic centre, y in the direction of 
Galactic rotation, and z towards the north Galactic pole, 
He and Hb denoting the proper motions in galactic longitude 
and latitude in mas, it the parallax in mas and k=4.7405 
the numerical factor that gives Vt in km -1 . The tangential 
velocity of a star in the plane of the sky is then defined by 



V T = 



sin £ cos bfi£— cos £ sin b fib 
cos £ cos b fie — sin £ sin b fib 
cos b fib 



(1) 



The tangential velocity can also be written in terms of 
the space velocity v (with e giving the direction to the star) 
as 



v - v R 



(I 



(2) 



Combining this we get the equation for the vertical ve- 
locity towards the north Galactic pole 



Kflb 



+ u cos I tan b + v sin I tan b 



(3) 



7TCOS& 

The ensemble of tracer stars then gives an estimate of 
the vertical velocity distribution function /(|wo|). Ideally, 
one would like to determine /(|u>o|) from the proper mo- 
tions of all the stars in the two samples. We found it more 
practical to determine /(|wo|) from stars at low Galactic 
latitude, since we use proper motions for the stars and ig- 
nore the (incomplete) radial velocity data for the sample 
and the velocity distribution function is dominated by the 
proper motions rather than the unknown radial velocity at 
low Galactic latitude. We chose to use a Galactic latitude 
cutoff of |&| < 12°. The choice of 12° galactic latitude was the 
best compromise between obtaining more stars by going to 
higher latitudes, while keeping |6| low in order to minimize 
the effect of the unknown radial velocity on /(|wo|)- 

There are several error sources that affect the velocity 
distribution function, /(|u>o|). An obvious one is that what 
we are measuring is not the distribution at the plane /(|wo|) 
but the mean of the distribution /(|w|) within the region 
delimited by |6| < 12°. Since samples taken further away 
from the plane are more dominated by kinematically hotter 
stars, this effect widens the velocity distribution. Another 
effect stems from measurement errors in the parallaxes — 
these scatter into the sample more distant stars. These stars 
must, in order to meet the apparent magnitude limit of the 
survey, be more luminous and hence younger and kinemat- 
ically colder, which would narrow the velocity distribution. 
Finally, we are only measuring the part of the vertical veloc- 
ity that comes from the proper motion in latitude. However, 
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Figure 3. Histograms of vertical velocity at the plane, /|uio|, for 
the A-star (dotted) and F-star (solid) samples. 

this is the smallest error of the three. From studies of the 
effect on /(|wo|) of varying |6| and the distance limit of the 
sample in simulated Hipparcos catalogues we estimate that 
the combined effect from these errors are of the order of 
0.1 km s -1 for the A-sample, and even smaller for the F- 
sample, at a cutoff latitude of \b\ < 12°. For |6| < 20° the 
correction rises to 0.5 km s -1 which would eventually lead to 
a 20 per cent error in the estimated local density of the disc. 
For |6| < 2° the uncertainty in /(|wo|) due to the declining 
sample size is also 0.5 kms -1 . Hence |6| < 12° represents 
the best compromise between declining sample size and the 
corrections described above. The measured velocity distri- 
butions for the two samples are also slightly broadened due 
to errors in the parallaxes and proper motions. However, 
simulations show that this effect is quite small for our sam- 
ples. The sample dispersion increases by a mere 0.04 km s -1 
for the F-stars, and by 0.16 km s -1 for the A-stars due to 
these errors. 

The velocity distributions for the two samples, /(|u>o|), 
are shown in Fig. 3, after a correction for broadening and 
for the Solar motion relative to the sample stars of u = 
10kms _1 ,i; = lOkrns -1 and w — 7kms~ 1 . This Solar mo- 
tion is based on an analysis of the the complete survey and 
is in good agreement with other determinations based on 
Hipparcos data (Dehnen & Binney 1998; Bienayme 1999). 
There are 723 A-stars and 683 F-stars in these two distri- 
bution functions. The velocity dispersion of the sample cor- 
rected for measurement errors and outliers is 5.7±0.2 km -1 
for the A-stars and 8.3±0.3 km -1 for the F-stars. 



3 DISC MASS MODEL 

We determine the local density of matter by comparing the 
density law of a stellar tracer, predicted from it's vertical 
velocity distribution at the plane, with the observed density 
law measured by Hipparcos. We follow the method of von 
Hoerncr (1960) in which the distribution of absolute vertical 



Figure 4. Resulting K z -force law for three mass models with the 
same total midplane density of po = 0.103 Mq pc~ 3 , of which 
0.093 Mq pc~ 3 is in the disc components and 0.01 Mq pc~ 3 is 
in the dark halo. The solid line shows the reference model from 
Table 1, the dashed line model is entirely made up of cold gas 
(<r z = 4 km s _1 ) and the dotted line of old stars (a z = 20kms~ 1 ). 



Table 1. Disc Mass Model 



i Description Pi(0) o~i Si Note 

M Q pc~ 3 km s _1 M pc~ 2 



1 


H 2 


0.021 


4.0 


3.0 


C 


2 


HI(1) 


0.016 


7.0 


4.0 


C 


3 


HI(2) 


0.012 


9.0 


4.0 


c 


4 


warm gas 


0.001 


40.0 


2.0 


c 


5 


giants 


0.0006 


17.0 


0.4 


H 


6 


M v < 2.5 


0.0031 


7.5 


0.9 


H 


7 


2.5 < M v < 3.0 


0.0015 


10.5 


0.6 


H 


8 


3.0 < M v < 4.0 


0.0020 


14.0 


1.1 


H 


9 


4.0 < M v < 5.0 


0.0024 


19.5 


2.0 


H 


10 


5.0 < M v < 8.0 


0.0074 


20.0 


6.5 


H 


11 


M v > 8.0 


0.014 


20.0 


12.3 


X 


12 


white dwarfs 


0.005 


20.0 


4.4 


L 


13 


brown dwarfs 


0.008 


20.0 


6.2 


L 


14 


stellar halo 


0.0001 


100.0 


0.6 


L 



Notes. C : component constrained by column density 

H : component constrained by local density using Hipparcos 

L : component constrained by local density using star counts 

X : component constrained by column density using HST star 

counts 



velocity of the tracer stars at z = 0, /(|wo|), may be 

used to construct the tracer's vertical density profile in the 
gravitational potential $ of the disc (Flynn & Fuchs 1994). 

Models for the vertical potential of the Galactic disc 
are constructed in the manner of Bahcall (1984a, b,c) in 
which the local visible components of the disc are modeled 
as isothermal components, specified by a local density and 
a velocity dispersion. We have made a number of refine- 
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ments to the mass models of Bahcall because of observa- 
tional advances since 1984, in particular for the upper main 
sequence, the giant branch, white dwarfs, red dwarfs and 
brown dwarfs. The basic model is shown in Table 1, and the 
resulting K z -force law is shown in Fig. 4 together with two 
simple one-component models with the same local density 
as the reference model. 



3.1 Upper main sequence and giant branch 

For the main sequence above the turnoff and for giants, Hip- 
parcos now provides much improved measurements of the 
stellar luminosity function. Holmberg et al. (1997) have de- 
veloped a version of the Bahcall-Soneria Galaxy model with 
revised luminosity /colour distributions and scale heights 
which fits the Hipparcos Survey data. We have determined 
the local number density of stars in the absolute magni- 
tude ranges used by Bahcall, for convenient comparison with 
his models. This affects rows 5 to 9 in Table 1. From our 
model fitting, there are 0.0013 stars pc" 3 for My < 2.5, 
0.0010 stars pc" 3 for 2.5 < M v < 3.0, 0.0015 stars pc" 3 for 
3.0 < M v < 4.0 and 0.0021 stars pc" 3 for 4.0 < M v < 5.0 . 
For the luminosity range 5.0 < M v < 8.0 we take 0.0090 
stars pc" 3 from Jahreifi & Wielen (1997). Adopting the 
mass-luminosity relation of Henry & McCarthy (1993), we 
derive the local mass densities pi shown in rows 6 to 10 
in Table 1. The local number density of giants is 0.0005 
stars pc . Adopting a mean mass of 0.9 we derive the lo- 
cal mass density of giants of 0.0006 M pc" 3 (row 5 of Ta- 
ble 1); Jahreifi & Wielen (1997) derive the same value also 
from Hipparcos. While the number densities of these com- 
ponents of the model can be accurately determined, there 
is a large relative error of order 20 per cent in the mass- 
luminosity conversion in determining individual pi. How- 
ever, since these stars contribute only about 10 per cent of 
the local mass, this error has only a minor impact on the to- 
tal mass model. For these stars the velocity dispersion can 
be estimated directly using the same technique as for the 
tracer stars (section 2.2) and are shown as the Oi in Table 1 
in rows 5 to 9. Note that the (Hipparcos based) observational 
data essentially constrain the local density and velocity dis- 
persion of these stellar types, and this is emphasised by the 
H in column 6 of Table 1. 



3.2 M dwarfs 

The column density of disc M dwarfs (Table 1, row 11) can 
now be measured directly via star counts using HST (Gould, 
Flynn & Bahcall 1998 and references therein) , rather than by 
extrapolation of the mass function. They measure a column 
density of M disc dwarfs of 12. 3± 1.8 Mq pc" 2 . We represent 
the M dwarfs as a single component with a velocity disper- 
sion of 20kms _1 , and we adjust the local density in each 
model we run to be consistent with a column density of M 
dwarfs of 12.3 M pc" 2 . These stars are better constrained 
observationally by column density than by local density, as 
indicated by an X in column 6 of Table 1. 



3.3 White dwarfs 

Several new estimates of white dwarf (WD) number den- 
sity have been made since 1984. Oswalt et al. (1996), in an 
analysis of white dwarfs discovered as proper motion com- 
panions to main sequence dwarfs, after allowing for for the 
fraction found in binaries, report the space number density 
7.6t 3 '7 x 10" 3 pc" 3 . For an adopted WD mass of 0.6 M 
we have 0.0046 M Q pc" 3 . Leggett, Ruiz & Bergeron (1998) 
find for single WDs in the proper motion survey of Liebert, 
Dahn & Monet (1988) a space density of 3.4 x 10" 3 pc" 3 
or 0.002 M pc" 3 . Jahreifi & Wielen (1997) report 0.005 
Mq -3 for 4 WDs confirmed by Hipparcos to be within 5 
pc, very close to the Wielen (1974) value of 0.007 M pc -3 , 
which can be compared with that adopted by Bahcall (1984), 
0.005 M pc" 3 . Knox, Hawkins & Hambly (1999) finds 
4.16 x 10" 3 pc" 3 or 0.0025 M Q pc -3 . All these estimates 
are lower limits, since they are based on identified white 
dwarfs. Interestingly Festin (1998) reports a much higher 
value of 0.013 M pc" 3 for 7 white dwarfs found using the 
innovative technique of searching against the opaque screens 
formed by dark molecular clouds in Orion and Ophiuchus. 
Considering the range of values reported, for the purpose of 
building mass models we adopt 0.005 M pc -3 but we will 
consider model WD densities which cover the range [0.002, 
0.013] M pc . WDs are primarily constrained by local 
density measurements as indicated by an L in column 6 of 
Table 1. 

3.4 Brown dwarfs 

The long-sought freely floating brown dwarfs (BD) of the old 
disc have now been found by four groups, namely the Calan- 
ESO proper motion survey (Ruiz, Leggett & Allard 1997), 
the DENIS mini survey, (Delfosse et al. 1997), a BRI survey 
of Irwin, McMahon & Hazard (1991) and UK Schmidt sur- 
vey of ESO/SERC field 287 (Hawkins et al. 1998). Fuchs, 
Jahreifi & Flynn (1998) summarize the findings of these 
groups and the implied minimum densities of BDs in the 
local disc. The four surveys lead to density estimates of BDs 
of 0.46 pc" 3 , 0.076 pc -3 , 0.069 pc" 3 , and 0.15 pc" 3 respec- 
tively. Adopting a BD mass of 0.065 M this gives mass 
densities of 0.03 M pc" 3 , 0.0049 M pc" 3 , 0.0045 M 
pc" 3 and 0.01 M pc" 3 , respectively. These surveys are all 
in their early stages and promise an accurate accounting of 
BDs within a few years. A fifth survey, 2MASS, has also re- 
ported at least 5 BDs in a 420 deg survey (Reid et al. 1998) 
although it is not possible yet to determine space densities 
for these BDs. In the model we adopt a BD mass density 
of 0.008 M pc -3 , and consider models over the range of 
[0.004, 0.03] M pc -3 . BDs are primarily constrained by lo- 
cal density measurements as indicated by an L in column 6 
of Table 1. 

3.5 Interstellar matter 

The remaining component in the model is by far the least 
well understood. Many studies and compilations exist in the 
literature on the structure and composition on the multi 
phase ISM (Hollenbach & Thronson 1987; Combes 1991) 
which describe a complex mixture, ranging from cold (10 
K) and dense (50 M pc" 3 ) dark molecular clouds, to the 
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Figure 5. The velocity distribution function for the stellar com- 
ponents of our model (solid line) compared to that observed in the 
solar neighbourhood by Fuchs & Jahreifi (private communication, 
histogram) . 

very hot (10 6 K) and dilute (0.0001 M pc~ 3 ) X-ray emit- 
ting coronal gas. Here we adopt the multi-phase model of 
Bahcall ct al. (1992) recognizing that the best determined 
parameters are the column density and velocity dispersion, 
which are observable quantities, whereas the deduced local 
density suffers considerably uncertainties (indicated by a C 
in column 6 of Table 1). This four-component model con- 
sists of 3 Mq pc~ 2 at 4 km s" 1 of molecular gas (Scoville 
& Sanders 1987), two neutral atomic components each of 
4 M pc~ 2 , one cold at 7 km s _1 and one warm at 9 km 
s _1 and finally one hot ionized component with 2 Mq pc~ 2 
at 40 km s -1 (Kulkarni & Heiles 1987). The total surface 
density of these components is about 13 Mq pc~ 3 , with an 
uncertainty of about 50 per cent. 

3.6 Adopted velocity dispersions 

For the early stellar types and giants in the model (rows 5 
to 9) the local velocity dispersion Ui is rather accurate be- 
cause they are derived from Hipparcos data. For all other 
non-gaseous components the velocity dispersion is not so 
well determined, and we simply adopt a uniform value of 20 
kms -1 for what amounts to K to M dwarfs, white dwarfs 
and brown dwarfs. We are concerned in this paper with us- 
ing the Hipparcos data to measure the local density of mat- 
ter, whereas the exact choice of Ui for these components 
only affects the column density of the models. To illustrate 
that the adopted value of 20 kms _1 is consistent with known 
constraints, we plot in Fig. 5 the velocity distribution of the 
combined stellar components in the model (rows 5 to 12 in 
Table 1) and an observational determination of /(|u>o|) by 
Fuchs & Jahreifi (private communication) based on Hippar- 
cos data for nearby stars. For purposes of determining the 
local density the match between model and data is suitably 
accurate. 



4 MODEL FITTING AND RESULTS 

We now use the velocity distribution functions, /(|ioo|), as 
determined in section 2.2 for the A- and F-star samples, to 
predict the density falloff in z of such stars in our disc mass 
model. The density falloff v(z) is determined via the relation 
(Fuchs & Wielen 1993; Flynn & Fuchs 1994): 

■•(-)-° r /(|wo|)wodw ° (4) 

j^rm \Jw\ - 2$ 

where wo is the vertical velocity at the mid-plane, z — 
and &(z) is the total gravitational potential generated by the 
mass model. Equation 4 ignores the radial term in Poisson's 
equation. If the term is included, Poisson's equation takes 
the form (Binney & Merrifield 1998): 

4-Gp=-^ + 2(B 2 -^ 2 ) (5) 

where A and B are the Oort constants. Using 
the values derived by Feast & Whitelock (1997), A = 
14.82 km s^kpc" 1 and B = -12.37kms _1 kpc" 1 the cor- 
rection term amounts to — 0.0025 Mq pc~ 3 . However, in 
a new determination of Mignard (1998) quoted in Ko- 
valevsky (1998), the values are A — ll^kms^ 1 kpc -1 and 
B = -lO.Skms^kpc -1 giving a smaller correction of 
— 0.001 Mq pc . These corrections are much smaller than 
other sources of error, and can be safely neglected. 

We integrate Eqn. 4 numerically using the model in 
Table 1 and the two velocity distribution functions /(|u>o|) 
shown in Fig. 3 for the A- and F-stars. The integration is 
performed by using a velocity interval of Awo = 0.01 kms -1 
in order to produce a fairly smooth number density falloff 
- a too large velocity interval (e.g. 0.1 kms -1 ) causes a 
visible and non-physical sawtooth effect in the the resulting 
number density law v. 

4.1 A-star sample 

The number density falloff for the A-stars resulting from 
the mass model in Table 1 is shown by the dotted line in 
Fig. 6, versus the Hipparcos data (histogram) and is already 
seen to be a good fit. We attempted to improve the fit by 
adding or removing mass from the model, and minimising 
the difference between the predicted and observed number 
density as a function of z using a standard \ 2 statistic. It 
made little difference to the fitting whether this mass was 
removed proportionally from all rows in the model or (more 
reasonably) from the most uncertain rows in the model, as 
one expects since the measurement is so local. For the A- 
stars the resulting best fitting local density was po = 0.103± 
0.006 Mq pc~ 3 . i.e. just the same as the reference model of 
Table 1. 

4.2 F-star sample 

We performed a similiar fitting for the F-stars as for the 
A-stars. For the F stars the best fitting vertical density dis- 
tribution is shown by the dashed line in Fig. 7, and corre- 
sponds to a local density of po = 0.094 ± O.O17M0pc -3 . 
This best fitting model has slightly less local matter (by ap- 
proximately 9 per cent) than the basic model in Table 1. As 
for the A-stars, it made little difference whether the mass 
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Figure 6. The vertical distribution of A-stars (histogram) ver- 
sus the prediction of the best fitting mass-model with po = 
0.103 Mq pc -3 (dotted line). Also shown is the F-star model with 
po = O.O94M pc~ 3 (dashed line). 
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Figure 7. The vertical distribution of F-stars (histogram) ver- 
sus the prediction of the best fitting mass-model with po = 
0.094 Mq pc~ 3 (dashed line). Also shown is the A-star model with 
po = 0.103M Q pc- 3 (dotted line) 

was removed in proportion from the whole model or from 
specific rows. To make this best fitting model mass was re- 
moved proportionally from the gas and white/brown dwarf 
components. 

4.3 Error estimates 

We have estimated our confidence limits via a series of 
Monte Carlo simulations of observations drawn from syn- 



thetic Hipparcos survey catalogues. The catalogues were cre- 
ated from a kinematical Galactic model constructed to fit 
the Hipparcos and Tycho catalogues (Holmberg 1999). In 
each simulation the number density and velocity distribu- 
tion of the tracer stars was exactly known. A sample was 
drawn from the simulation taking into account the Hippar- 
cos magnitude limits and magnitude dependent parallax and 
proper motion errors. The local density of the disc was then 
determined from the artificial sample in the same manner 
as for the actual data, by measuring v{£) and /(|wo|) and 
fitting mass models. Several thousand Monte-Carlo simula- 
tions of this type allowed us to estimate our confidence limits 
for a given disc density. 

In the modeling of the A-star sample the velocity dis- 
persion was 5.72 ± 0.18 km s -1 and the best fit density is 
po=0.101 ± 0.006 M Q pc~ 3 , with a 95 per cent confidence 
limit of ± 0.011 Mq pc~ 3 , showing that the uncertainty in 
the density is totally dominated by the quadratic depen- 
dence of the velocity dispersion. The F-star sample had a 
velocity dispersion of 8.14 ± 0.24 kms -1 with a best fit den- 
sity of po = 0.103 ± 0.017 M Q pc~ 3 , with a 95 per cent 
confidence limit of ± 0.023 Mq pc -3 . This shows a differing 
effect in which the variance of the resulting density is domi- 
nated by the uncertainty in the determination of the fall-off 
of the tracer density. 

In the x 2 fitting process the density distribution result- 
ing from the mass model is considered to be a distribution 
function without errors. This has the effect that the result- 
ing value of x 2 is increased by the errors from the velocity 
distribution function of the tracer stars. In the simulation of 
the complete sample the mean value of x 2 was 59.7 instead 
of the expected 56.3 for a model without errors. For the ac- 
tual Hipparcos samples the resulting values of x 2 was 77.76 
for 57 DF in the total sample, 44.59 for 31 DF in the A-star 
sample and 33.11 for 25 DF in the F-star sample. 

4.4 Disc column density 

The total column of gravitating matter in the mass models 
that best fits all our samples is quite close to the 48 Mq pc~ 2 
of visible matter in our basic disc mass model of Table 1. We 
strongly emphasize here that our samples are so local that 
we cannot put any constraints on the column, since most of 
it is above our tracers. An illustration of this is that if the 
velocity dispersion of the hot disc stars in the mass model 
is increased from 20kms~ 1 to 25kms~ 1 the column inflates 
by 7Mqpc~ 2 but the local density estimate changes with 
only 0.002 M pc -3 (see also section 3.6 above). 

4.5 Dependence on disc mass model 

How critical is the formulation of the mass model to the 
result? We tested this by using the quite different model 
from Flynn & Fuchs (1994), slightly refined by the inclusion 
of Hipparcos observations (Fuchs 1998, private communi- 
cation). This model is quite different in that rather than 
modeling individual mass components by isothermals, as we 
do here, Flynn & Fuchs reconstructed the disc potential $ 
using the velocity distribution function of K and M dwarfs, 
to which a small component of young stars and gas had 
been added. Adopting this model for $ resulted in a change 
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Table 2. Local density estimates 



Density po. 




Reference 


(M pc- 3 ) 


(M 0P c- 3 ) 




0.185 


0.020 


Bahcall (1984b) 


0.210 


0.090 


Bahcall (1984c) 


0.105 


0.015 


Bienayme, Robin & Creze (1987) 


0.260 


0.150 


Bahcall et al. (1992) 


0.110 


0.010 


Pham (1997) 


0.076 


0.015 


Creze et al. (1998) 


0.102 


0.010 


This paper 


0.150 


0.026 


Straight average 


0.108 


0.011 


Variance weighted average 



of only 0.003 Mq pc 3 to our best-fitting local mass density 
determinations. 



4.6 Combined A- and F-star sample 

Our conclusion in this section is that the dynamically es- 
timated mass in the solar neighbourhood is completely ac- 
counted for by the identified material in gas, stars and stellar 
remnants without any need for dark matter in the disc. Since 
both the A- and F-star samples agree on the local density 
within their errors, we combined the sample and derived a 
best fitting local density of p =0.102 ± 0.006 M pc~ 3 , with 
95 per cent confidence limits of ± 0.010 M pc -3 . 

4.7 Comparison to other work 

Our estimate of the local mass density is po = 0.102 ± 
0.006 M pc~ 3 . From very similar Hipparcos data Pham 
(1997) estimates p = 0.11 ± 0.01 M pc~ 3 and Creze et 
al. (1998) estimate p = 0.076 ± 0.015 M pc" 3 . These esti- 
mates differ at less than the 2— a level, but we have investi- 
gated possible causes for the discrepancy between them. The 
major difference between the determinations is the method 
applied to the data. Pham (1997) used the 2nGh 2 z po = a z 
relation, which implicitly assumes that the tracer stars (i.e. 
the A- and F-stars) have a Gaussian velocity distribution. 
Fig. 3 clearly illustrates that this is in fact only an approx- 
imation. Had we carried out our analysis using a Gaussian 
distribution function rather than the measured distribution 
function /(|u>o|) we would have recovered a local density of 
po = 0.110 M pc -3 . Creze et al. (1998) estimate the local 
density by representing the vertical acceleration in the disc 
K z by a linear function, K z = az, where a is a constant. K z 
may not actually be linear in the region of interest (z < 100 
pc). As can be seen in Fig. 4., K z would not be well rep- 
resented by a linear function in the model adopted here for 
all known mass components of the disc (solid line). If we 
had assumed a strictly linear vertical force law, then tests 
showed that our best fitting local density would have been 
reduced to 0.089 M pc~ 3 . This effect could explain part of 
the difference between our result and Creze et al. (1998). In 
summary, the small differences between the three local den- 
sity estimates using Hipparcos data are probably a result of 
the three different methods used. 



Stothers (1998) gives a detailed summary of measure- 
ments of the local disc density over more than 60 years since 
Oort (1932). Stothers draws attention to the fact that there 
are sufficiently many determinations of po that one might 
apply the central limit theorem and derive the best density 
from the straight average, which is 0.141 ± .0.010 M pc~ 3 . 
However, in contrast to most older studies, almost all deter- 
minations from the last 20 years have clearly stated uncer- 
tainty estimates attached to the preferred density. The data 
can be found in Table 2 together with our own estimate in 
this paper for po- The table shows that high estimates of the 
density generally have large uncertainties. The sample mean 
weighted by the inverse variance is 0.108 ± .0.011 M pc~ 3 
which is quite consistent with the Hipparcos results. One 
should note that the last three determinations are based on 
almost the same data from Hipparcos. 



5 CONCLUSIONS 

From a sample of A and F stars taken from the Hipparcos 
catalogue, we determined vertical space velocities and den- 
sity distributions. Using a method of von Hocrner (1960) we 
solved for the vertical potential and hence local density that 
would give the observed velocity and density distributions. 
We find that local dynamical mass density of the solar neigh- 
bourhood is 0.102 ±0.010 M pc -3 well compatible with the 
identified material in ordinary matter leaving no need for 
dark matter in the disc. Our estimate is found to be in rea- 
sonable agreement with other estimates, taking errors into 
account. We conclude that there is no compelling evidence 
for significant amounts of dark matter in the disc. 
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